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1 .  Introduction 

Recently,  the  temperature  variations  of  various  sound  propagations  have 

3  i  ^  2  3  ^ 

been  an  important  subject  in  the  studies  of  liquid  "He,  4He,  JHe-  He 

mixtures^  and  solutions  of  ^He  in  liquid  ^He ,  ^  because  their  properties  are 

3 

closely  related  to  the  exitation  spectrum  of  liquid  He  and  he,  which  can 
determine  the  transport  coefficients^  through  the  kinetic  equations.  It  is 
well  known  that  at  low  temperatures  and  low  pressures  the  thermal  properties 
of  superfluid  helium  are  not  determined  mainly  by  the  rotons  but  rather  by  the 
low-momentum  acoustic  phonons.  These  phonons  do  not  possess  the  normal 
dispersion  relation,  but  instead  exhibit  an  anomalous  excitation  spectrum.  In 
fact,  these  two  cases  are  quite  different  in  microscopic  processes.  The 
former  contributes  to  four-phonon  processes  (4PP)  while  the  latter  is  governed 
by  three-phonon  processes  (3PP) . 

Recently,  Andreev  and  Khalatnikov  have  evaluated  the  temperature 
variation  of  first  sound,  and  Singh  and  Prakash7  have  used  the  retarded 
single-particle  Green's  function  for  a  weakly- interacting  Bose  gas  to  obtain 
the  first  sound  by  using  the  wrong  normal  dispersion  relation.  However,  in  a 

g 

recent  paper  we  have  developed  the  Landau-type  elementary  excitation 
spectrum,  which  is  anomalous  phonon- like  at  low  momenta  and  roton-like  at  high 
momenta  in  two-  and  three-dimensional  liquid  helium.  Starting  with  this 
elementary  excitation,  we  have  successfully  derived  and  explained  not  only  the 

various  sounds^  and  sound  attenuations,^  but  also  thermal  conductivity^  and 

12  U 

viscosity  in  partly  one-,  two-  and  three-dimensional  liquid  He. 

Regarding  the  evaluation  of  first  and  second  sound, ^  we  have  used  a  new 

approach,  which  takes  into  account  a  collision  term  in  the  Boltzmann  equation, 

and  have  obtained  first  and  second  sound  simultaneously.  We  reported  that  the 

second  sound  in  thin  helium  films  is  approximately  2  times  the  first  sound 


3 


velocity  in  the  low- temperature  and  low-frequency  limit.  At  low  frequencies 

such  that  u>  r  <<  1,  where  u>  is  the  sound  frequency  and  r  is  the 

characteristic  time,  we  can  make  use  of  a  hydrodynamical  approach  to  sound 

propagation.  However,  for  the  opposite  case  of  w  r  »  1,  it  is  better  to  use 

s 

the  kinetic  equations.  In  this  collisionless  region  the  first  and  second 

sound  obtained  in  superfluid  hydrodynamic  equations  are  involved  in  the 

attenuation  coefficients,  which  contain  the  four  coefficients  of  second 

viscosity,  i.e.,  f ^ f  ^  and  f  These  four  coefficients  play  a  very 

important  role  in  the  investigation  of  sound  attenuation  in  bulk  liquid 

helium.  In  the  bulk  case,  the  attenuation  of  second  sound  depends  on  all  four 

viscosity  coefficients  together.  However,  the  term  which  generally  contains 

the  thermal  conductivity  strongly  affects  the  attenuation. 

Recently,  we  have  adopted  the  Landau  and  Khalatnikov  theory  to  derive 
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the  thermal  conductivity  and  viscosity  for  three  ranges  of  temperature, 
T  <  0.3  K,  0,3  K  <  T  <  0.8  K  and  T  >  0.8  K,  in  which  the  scattering  depends  on 
the  nature  of  interactions  between  elementary  excitations.  For  T  <  0.3  K  it 
is  important  to  note  that  the  3PP  do  not  affect  the  thermal  conductivity  but 
cause  the  viscosity  to  have  a  T  ^-dependence,  which  was  not  proven  by  Landau 
and  Khalatnikov' s  results.  The  4PP  and  phonon- roton  scattering  govern  the 
whole  mechanism  in  the  range  of  temperature  0.3  K  <  T  <  0.8  K,  and  for  T  > 
0.8  K  the  5PP  and  the  phonon-roton  scattering  contribute  mostly  to  the 
transport  coefficients. 

The  main  purpose  of  the  present  paper  is  to  evaluate  the  four 

coefficients  f  ^  f  ^ an<*  of  the  second  viscosity  as  a  function  of 

temperature  by  solving  the  superfluid  hydrodynamics  for  the  above  temperature 

14 

ranges  through  the  theory  of  the  kinetic  phenomena  developed  by  Khalatnikov, 
and  especially  for  T  <  0.3  K  we  will  investigate  the  temperature  variation  of 
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Che  four  viscosity  coefficients.  Throughout  this  paper  we  will  use  the  two- 
dimensional  elementary  excitation  spectrum  obtained  by  the  microscopic  ring- 
diagram  approximation,^ 


2  4 

£(p)  -  cQp[l  +  7p  -  S p  + 


(1.1) 


£(P)  -  A  + 


<p-?o>2 

2"0 


(1.2) 


where  p  and  P  are  phonon  and  roton  momenta,  respectively,  c^  is  the  sound 
velocity  at  absolute  zero  temperature,  7  and  S  are  positive  constants  which 
can  be  determined  by  the  potential  parameters,  and  A,  Pq  and  Mq  are  the  roton 
parameters.  Here,  we  have  adopted  a  soft  potential  with  a  Lennard- Jones- type 
tail  which  helps  to  make  a  smooth  connection  between  the  attractive  part  and  a 
soft  repulsive  core: 


<H  r) 


r*  r/d\12  .3,6, 

W  'V  ] 


r  <  a 


r  >  a 


(1.3) 


In  this  paper  we  define  the  liquid  helium  film  as  two-dimensional  -- 
less  than  three  atomic  layers,  namely  one  statistical  layer  of  3.6  A  --  and 
neglect  substrate  effects.  In  Sec.  2,  to  obtain  the  temperature  variation  of 
Che  viscosity  coefficients  in  various  processes,  we  will  investigate  the 
absorption  and  emission  processes  between  elementary  excitations  through  the 
collision  integral  equations,  and  we  will  solve  the  super  fluid  hydrodynamic 
equations  to  obtain  four  coefficients  of  the  second  viscosity  in  Sec.  3. 
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Finally,  in  Sec.  4  we  will  give  the  results  and  discussion  together  with  the 
thermodynamic  properties  in  terms  of  a  table  and  graphs. 

2.  Absorption  and  emission  of  phonons  and  rotons . 

To  describe  the  interactions  between  elementary  excitations  which  cause 
liquid  helium  to  make  a  transition  from  one  state  to  another,  we  can  consider 
generally  three  collision  processes:  (a)  emission  or  absorption  of  phonons, 
(b)  transformation  of  phonons  into  rotons  and  vice  versa,  and  (c)  emission  or 
absorption  of  rotons.  In  addition,  we  have  the  3PP,  4PP  and  5PP.  In  the  case 
of  (c)  the  energy  of  one  roton  would  have  to  be  at  least  2A,  such  that  it  can 
decay  into  two  rotons  each  with  energy  of  about  A.  However,  this  kind  of 
three-roton  process6  is  highly  improbable,  so  that  we  may  neglect  it. 
Therefore,  from  now  on,  we  only  take  into  account  the  3PP  and  5PP  and 
transformations  of  phonons  into  rotons  and  vice  versa. 

We  consider  first  the  3PP  (p^  +  p  «-  p^).  Since  the  total  number  of 
phonons  traveling  in  a  given  direction  is  changed  by  the  small -angle  3PP ,  the 
distribution  function,  which  depends  on  the  chemical  potential  a'  ,  can  be 
written  as 

n-  [exp(a'+pc)/kBT  -  l]'1  .  (2.2) 

The  difference  between  Eq .  (2.2)  and  the  equilibrium  distribution  function  n^ 
can  be  expanded  in  a  power  series  to  give 


•  no  "  •  no(no+I)  a' 


<5n  -  n 


(2.3) 
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The  collision  integral  J^pp(n)  and  the  kinetic  coefficient  r3pp  are  related  to 
each  other  through 


k«T  I 


c  T  I  d3PP^  2  ”  aP3pp 

B  J  (2jt  h) 


(2.4) 


and  the  collision  integral  J^pp(n)  obtained  in  the  second  quantization 
method^  is  given  by 


J3PP<n> 


*c0  2  r 

W.  (u+1)  P1P5P  6n(nn  ~ nn  )  S  (E  f-E  ) 

^ p0  J  1  2  p10  P20  f  1 


,  (2.5) 


(27rh)' 


where  u  is  Griineisen's  constant,  and  n  represents  the  equilibrium 

P  iO 

distribution  function  for  the  phonons  with  momentum  p„  Performing  the 
integration  for  Eq.  (2.5)  over  momentum  space,  we  obtain 


J3Pp(n>pdp  “ 


8  *V3pr 


••  © 


(2.6) 


Comparing  Eqs.  (2.4)  and  (2.6),  we  can  easily  obtain  the  kinetic  coefficient 
f3pp  for  the  3PP  as 


ic  2  w5  6 
16*  p0*  c0 


<kBT)‘ 


(2.7) 


As  for  the  4PP  case,  the  direction  of  momenta  of  the  colliding  particles 
is  not  changed  in  the  5PP  case.  The  calculation  of  the  transition  matrix 
elements  by  second-order  perturbation  theory  is  very  complicated,  and  some 
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terms  contain  vanishing  denominators  under  the  condition  y  -  0  in  Eq.  (1.1) 
and  collinear  scattering.  The  5PP  has  maximum  probablitv  in  small-angle 
scattering  and  leads  to  equilibrium  for  the  phonons  in  the  direction  of 
motion.  Instead  of  these  tedious  calculations  using  second-order  perturbation 
theory,  we  make  use  of  the  kinetic  coefficient  given  by  Landau  and 
Khalatnikov^  to  obtain  the  f^pp  as^ 


r 

5PP 


aT 


(2.8) 


Here,  a  is  constant  which  can  be  determined  experimentally  by  the  attenuation 
coefficient  of  ultrasonic  waves. 

Now  we  consider  the  transformations  of  phonons  into  rotons  and  vice 
versa.  When  energetic  phonons  with  energy  on  the  order  of  A  collide  with 
rotons  (p^  +  P^  «-  P^  +  **4)’  w^ere  P  and  p  represent  the  phonon  and  roton 
momenta,  respectively,  the  transformation  can  occur.  Therefore,  a  phonon  has 
at  least  a  very  large  energy  A,  and  the  interactions  between  phonons  and 
rotons  are  very  similar  to  that  of  the  scattering  between  rotons.  Adopting 
Landau  and  Khalatnikov' s  assumption,  ^  we  may  take  the  interactions  between  a 
roton  and  phonon  to  be  a  5 -function  potential, 


V  -  VQ6(r1-r2)  ,  (2.9) 

where  is  an  interaction  constant,  and  r^  and  r^  are  the  position  vectors  or 
the  phonons  and  rotons,  respectively. 

The  rate  of  change  per  unit  time  in  the  roton  number  due  to  the 


transformation  of  rotons  into  phonons  is  given  by 
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Nr  -  -  JJI  {niN2(N3+l)(N4+l)  -  (n1+l)(N2+l)N3N4)dw 


dp1dP2 

(2^)4 


(2.10) 


“t 

where  n.  and  N.  are  the  distribution  functions  with  momenta  p.  and  P. , 
11  *ii 

respectively.  du  is  the  differential  decay  rate  for  the  transformation  of 

—♦  —¥  -4 

phonons  and  rotons  with  momenta  p^  and  P2  into  two  rotons  with  momenta  P3  and 
P^  and  is  given  by 


o  o  dP « dP . 

du  -  —  |v  |Z  S(E  +  E  -  E  -  E  )  — - — - 

W  'IF1  K' 1  2  3  4;  4  ’ 

*  (2*X) 


(2.11) 


where  Iv^l  is  the  transition  matrix  element,  and  E ^  and  E ^  are  the  phonon  and 
roton  energies,  respectively.  Combining  Eq.  (2.10)  with  (2.11),  we  get 


Nr  -  -  JJI  (n1N2(N3+l)(N4+l)  -  (n1+l)(N2+l)N3Ni 


2jl  I  t,  1 2  dpldP2dP3dP4 

X  V  I Vipl  HEl  +  E2  *  E3  '  E4  (2^)8 


(2.12) 


Now  we  construct  the  symmetrized  pairwise  plane  wave  over  incoming  and 
outgoing  phonons  and  rotons  as 


p  v  i_  r  i(prVf2,?2)  i(vVVV 

^  pi ’ P2  “  yf  (e  +  e 


.  1_  ,  J<(P3’rl+P4’r2)  M(^3’r2+P4*rl). 

0(p3.p4)  "  yf  (e  +  ef 


(2.13) 
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and  combining  Eq.  (2.13)  wich  (2.11),  we  can  obtain  the  transition  matrix 
elements  as 


V 


IF 


2V 


>n'*I 


ji<P3+Vpl'P2)-rl 


(2.14) 


Here,  dQ  represents  the  area  element.  Making  use  of  the  <5  -function  identity 
together  with  Eq.  (2.14),  we  obtain 


5(p3+p4-prp2) 


(2*K)‘ 


I 


^(p3+p4-P1-P2) *ri 


(2.15) 


and  performing  the  integration  over  momentum  space,  we  get 


dP„ 


(2 


(2.16) 


As  mentioned  earlier,  from  the  fact  that  the  roton  momentum  F  is  almost  the 
same  as  that  of  the  roton  parameter  Fq,  the  expansion  of  the  roton 
distribution  function  N  as  a  function  of  chemical  potential  yields  Eq .  (2.12) 
as 


N 


Mr  ~  Moh  16*2 


kBT 


X 


(2»X) 


^  II  *3«* 


40dP3dP4 


(2.17) 


where  and  /i  ^  are  the  roton  and  phonon  chemical  potentials,  respectively, 
and  performing  the  integration  over  momentum  space  we  arrive  at 
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N 

r 


V 


^|vQ|  V 

3  2 

*  c0 


(2.18) 


where  is  the  number  density  of  roton  given  by 


(2.19) 


The  rate  of  change  per  unit  time  in  the  roton  number  can  be  expressed  in  terms 

of  the  kinetic  coefficient  T  ,  as 

ph-r 


rph-r^r  ^ph') 


(2.20) 


Comparing  Eqs.  (2.20)  and  (2.18),  we  can  deduce  the  kinetic  coefficient  ^ 

for  the  transformation  of  rotons  into  phonons  and  vice  versa  as 


r 


ph-r 


**|v0|  V 

3  2 


(2.21) 


Since  is  not  known  experimentally,  we  can  simply  rewrite  Eq.  (2.21)  as 


ph-r 


b  e 


-2A/kBT 


(2.22) 


where  b  is  a  constant  which  contains  several  parameters  of  the  elementary 
excitation  and  can  be  determined  experimentally  by  the  ultrasonic  attenuation. 
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3.  Coefficients  of  two-dimensional  second  viscosity 

When  liquid  ^He  is  in  its  equilibrium  state,  the  distribution  functions 
for  phonons  and  rotons  are  expressed  by 


n0ph 


>xpC 


p,(Vvn}  -1 

7t  1  -  13 

V 


(3.1) 


Or 


'  exp  { 


E  -p*(v  -v  ) 

- H-S-j,*! 

kBT  1J 


(3.2) 


However,  due  to  the  presence  of  the  energy  dissipation,  nonequilibrium  effects 
appear  in  energy  and  momentum  conservations,  superfluid  flow,  energy  flux,  and 
especially  entropy  which  is  not  conserved  but  increases  (entropy  increases  in 
order  to  determine  the  unknown  dissipation  coefficients).  Taking  into  account 
all  of  these  conditions,  we  can  accumulate  all  these  effects  within  the 
following  two  hydrodynamic  equations  which  concern 


SI 


£  +  VP  +  v^Vj  +  (j  -  pv^)Vv^  +  (j*^)vs  +  (vn,,;7)  ( j  'Pvs) 


V{f1'7(j -pvm)  +  C2Vvn) 


(3.3) 


<3v  v 

jf  +  +  i)  ’  + 


(3.4) 


Here,  j  and  p  are  the  momentum  and  mass  density,  P  is  the  pressure,  /i  is  the 
chemical  potential  of  the  liquid  helium,  and  C  ^  anc*  are 
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coefficients  of  second  viscosity.  The  coefficients  and  f  are  equal 
according  to  the  Onsager's  reciprocity  principle. 

It  is  obvious  that  the  coefficients  of  the  second  viscosity  depend  on 
the  variations  of  the  temperature  and  thermodynamic  functions,  and  thus  the 
number  of  phonons  and  rotons  are  changed  by  the  various  processes  that  we  have 
mentioned  earlier.  Let  N  ^  and  N^_  be  the  number  of  phonons  and  rotons  per 
unit  area,  respectively,  and  /i  ^  and  Mr  be  their  chemical  potentials.  When 
the  system  deviates  slightly  from  its  equilibrium  state,  and  N  ^  are 
changed  in  time  and  try  to  return  to  their  values  in  the  equilibrium  state. 
We  assume  that  the  nonequilibrium  distribution  function  n  deviates  very 
slightly  from  equilibrium,  i.e.,  small  deviations  of  the  density  and  entropy 
can  be  determined  by  the  time  derivatives  of  the  phonon  and  roton  numbers,  N  ^ 
and  Nr,  in  the  expansion  of  the  chemical  potential.  Since  the  phonons  and 
rotons  take  part  in  normal  motion  with  velocity  v^,  neglecting  the  quadratic 
effects  and  taking  only  linear  terms  in  and  we  can  obtain  the 
following  equations: 


N  +  N  Vv  —  -  7  u.  +  ~t  ,  u 
r  r  n  rr  r  'ph-r 


ph-r  ph 


(3.5) 


N  ,  +  N  ,Vv  —  7  ,  -  7  n  , 

ph  ph  n  ph-r  pp  ph 


(3.6) 


The  kinetic  coefficients  on  the  right-hand  sides  of  Eq .  (3.5)-(3.6)  are 

symmetric  in  the  indices  r  and  ph. 

Since  the  rate  of  change  per  unit  time  in  the  density  p  and  entropy  S 
can  be  expressed  by  the  continuity  equations 


p  +  v.j  -  0 


(3.7) 
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S  +  VSv  =>  0 
m 


(3.8) 


substitution  of  Eqs .  (3.7)-(3.8)  into  (3.5)-(3.6)  yields 


N 

r 


i'7pv 

n 


oN 

~ — “  p7v 
op  n 


3N 


2h 


oS 


SVv  + 
n 


oN  , 

— Vpv 
op  n 


oN  , 

op  n 


(3.9) 


(3.10) 


Combining  Eqs.  (3.  9) -(3.10)  with  Eqs.  (3. 5) -(3. 6),  we  can  express  and  N  ^ 
in  terms  of  p  and  S  bv 


3N  oN  3N  _ 

3?  +  (Nr  •  '  ?TP)VVn  ’  ‘7rrMr  +  7ph-rMph  1  (3>11) 


3N  ,  3N  3N 

~fi~  +  (Nph  •  ^"S  -  ^)7vn  ‘  7ph-rV7ppMph 


(3.12) 


The  kinetic  coefficients  7  ,  7  ,  and  7  15,17  m  be  replaced  as 

rr  ph-r  pp  -  r 

r-  ,  T_  ,  r  ,  and  ro  ,  which  are  the  kinetic  coefficients  of  the  three- 
3rp  5pp  ph-r  3pp 

roton  process  (3RP),  5PP,  phonon-roton  interaction  and  3PP,  respectively,  as 
follows 


7  ■  r.  +  r  , 

rr  3rp  ph-r 

*y  *  y  »  r 

r-ph  ph-r  ph-r 
7pp  "  r3pp  +  r5PP  +  rph-r 


(3.13) 
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Since  is  very  small  compared  to  the  other  kinetic  coefficients  ,  we  may 

neglect  it.  Substituting  Eq.  (3.13)  into  Eqs  .  ( 3 . 11)  -  ( 3 . 12)  and  solving  for 

M  ,  and  n  ,  we  obtain 
ph  r 


"ph  "  r3PP  *  r5PP  '  <N  ■  fiS  •  f, ’f)  VVn' 


(3.14) 


*r  -  ? - (f  5(J-^n)  ■  (K  -  ffs  .  fp)  V? •  , 

r  1  3pp  +  ^pp  °p  n  as  ap  n 


,  3N  3N  3N 

+  {~  ^ '  ^p>  Vv"’  • 


(3.15) 


where  N  -  N  +  N  ,  . 

r  ph 

In  Eqs.  (3. 3) -(3. 4)  the  pressure  P  and  chemical  potential  m  can  be 
expressed  in  terms  of  the  chemical  potentials  j*  ^  and  as 


(3.16) 


+  jj£>> 


Solving  Eqs.  (3 . 14) - (3 . 16)  together  with  Eqs.  (3.3)-(3.4)  for  the  coefficients 
f 1 .  f 2 >  ^3  an<3  of  the  second  viscosity,  we  obtain 


c  d)  -  -  1  ap  aN  1  ap  r 

1  r3PP  +  r5PP  dfir  aV  dp  '  rph-r  3"r  3p 


(3.17) 


r  m  _  - 1 _  ,32-  .  i£_  Me  M  .  1  c?_ 

r2(T)  +  rv  3u  +  au  .  )(N  *  ass  ‘  a„p)  +  r  .  au 


3PP  +  r5PP  dfir  df*ph 


as  dpr/  r  ,  dn 

ph-r  r 
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<3N  3N 

x  (N  +  zr^s  -  r-^p) 
r  oS  ao 


(3.18) 


M(t)  “  '  r - T~r -  + 

3  r3PP  +  r5PP  ^r  *Mph  dt> 


i  -  °N 
1  au  _ r 

r  ,  <5m  dp 

ph-r  r 


(3.19) 


f,(T)  -  - ~r —  (f^-  +  z^-)(N  -  ?fs  -  f^p)  +  z-^ — 

4  r3PP  +  r5PP  3**r  3**ph  oS  3p  rph-r  apr 


3N  3N 

*  (Nr  *  ifS  *  jfp) 


(2.20) 


Making  use  of  the  thermodynamic  identities 


dE_  -  TdS  +  pd p  -  N  d^t  -  N  ,  d^t  ,  , 

0  r  'r  ph  ph 


(3.21) 


P  -  -Eq  +  ST  +  pp  , 


(3.22) 


we  can  transform  Eqs.  (3 . 17) - ( 3 . 20)  into  the  following  forms: 


i  a»  5N  3N 

ri(T)  “  r3pp  +  r5pp  a7(N  '  ass  ‘  ~Pp) 


.  3N  SN  3N 

zr^ —  z- ^(N  -  ~s  -  r -f-p) 

rph-r  ap  r  °S 


(3.23) 


.  3N  3N  .  3N  3N 

r,(D  -  Z - TT -  (N  *  -  —  p)1  +  (N  -  ZZ^S  -  jfp)  ,  (3.24) 

2  r3PP  +  r5PP  5S  3p  Fph- r  r  °S  °P 


r  (T)  -  - - -  [— 1  ^  1 — 

3  r3PP  *  r5PP  W  rph-r 


.  /-3N  \  , 

bh; 


(3.25) 
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f4(T)  - 


r  +  r 

3PP  5PP 


3N, 


3N  3N 


?;(N  ■  iis  - 


ph-  r 


5N 

_ r 

dp 


3N  a  N 

X  (Nr  *  ifS  *  jfp)  -  ^1(T) 


(3.26) 


Here,  and  C/  are  shown  co  be  equal  as  expected  from  Onsager's  reciprocity 

symmetry  relation  for  the  kinetic  coefficients. 

As  we  have  recently  proven,  near  T  -  0  K  the  3PP  are  the  main  influence 

11  12 

on  the  first  viscosity  coefficient  in  2D  and  3D  liquid  helium,  which  was 

•  10  . 

not  shown  by  Landau  and  Khalatnikov,  and  in  ultrasonic  sound  attenuation  tne 
3PP  also  influence  the  coefficients  of  second  viscosity  near  T  -  0  K. 
Therefore,  to  investigate  the  contribution  of  the  3PP  to  the  coefficients  of 
the  second  viscosity,  we  should  consider  the  3PP  effects  separately  from  other 
processes.  To  do  this,  we  only  keep  the  terms  for  the  3PP  in  Eqs.  (3.3)* 
(3.12)  and  then  obtain 


3N 


ph 


3PP 


dp 


*  •  "V  •  <\h 


3N  ] 

as 


iS  - 


3N 

y~p )  Vv  i 
dp  n 


(3.27) 


The  pressure  P  and  chemical  potential  4  in  Eqs.  (3.3)-(3.4)  and  Eq.  (3.6), 
which  depend  only  on  the  phonon  chemical  potential,  can  be  written  as 

VP  -  „  1  V4  -  n  1  (3.28) 

i^ph  phJ  Kh  w 

Therefore,  with  the  use  of  the  thermodynamic  identities  Eqs.  ( 3 . 21)  -  ( 3 . 22)  ,  we 
can  obtain  the  contributions  from  the  3PP  to  the  coefficients  of  the  second 


viscosity  in  Eq .  (  3 . 2 3 ) - (  3 . 26 )  ,  together  with  Eqs.  (  3 . 2 7  )  -  (  3  .  2  8  )  ,  near 

T  *  0  K  as 


fx(T) 


r2(T)  - 


r3(T) 


f4(T> 


1 

3Nph 

r 

3PP 

dp 

l 

r 

3PP 

>v  • 

l 

r!!W 

r 

3PP 

[  So  J 

_ l_ 

<3N  , 
ph 

r 

3PP 

dp 

'  are 

4 

equal 

5N  3N  , 

_ Eh  _ 2h  ■ 

N  u  -  «c  S  ■  „  P 

ph  oS  op 


a N  ,  3N  , 

_ Eh<-  _ Eh  i  2 

3S  S  dp  p]  ’ 


(3.29) 


(3.30) 


(3.31) 


3N  .  3N  , 

(NPh  *  ' -t-* 


(3.32) 


principle . 


4.  Results  and  discussion 

In  the  previous  sections  we  have  derived  the  kinetic  coefficients  for 

two  temperature  regions,  and  making  use  of  these  coefficients  we  have  obtained 

the  four  coefficients  of  the  second  viscosity.  For  T  <  0.3  K  the 

contributions  of  the  3PP  to  the  coefficients  are  most  important,  and  for  T  > 

0.3  K  we  may  neglect  the  contribution  from  the  3PP.  Therefore,  dropping  the 

r3pp  terms  in  Eqs.  ( 3 . 23)  -  ( 3 . 26)  ,  these  coefficients  of  the  two-dimensional 

16 

(2D)  second  viscosity  are  reduced  to  those  of  the  bulk  case,  except  tor 
differences  between  2D  and  3D  thermodynamic  dimensionality.  Therefore,  we 
will  not  write  the  coefficients  here  repeatedly  for  T  >  0.3  K. 
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To  analyze  concreteLy  the  four  coefficiencs  of  the  2D  second  viscosity 

■  ,  .  8 

for  T  <  0.3  K,  we  use  the  following  two-dimensional  thermodynamic  functions: 


kR  Trm  9  5 !  r  ( 5)  B 

VT>  -  T.  <kBT)  +  — —  <v>  * 


(4.1) 


Sr<T,  -  p^j 


M0kBTl H  kBP0  ,3 


,3  A  -VV 

2  +  k  T  e 
B 


(4.2) 


vj  m  _  L121  fV]2 
Ph  2jt  [pc  J 


(4.3) 


«  *•  <j 

Here,  B^  -  [1  +  ~  (~a- v^) ] /2Cg ,  n  is  the  number  density  of  2D  liquid 

helium,  and  the  other  parameters  are  given  in  Eq  .  (1.3).  To  see  the 

temperature  variation  of  the  coefficients  numerically,  we  adopt  the  parameters 

1 8 

which  are  deduced  from  the  specific  data  of  Bretz  et  al .  To  explain  the 

9  10 

various  experimental  data,  ’  we  have  made  use  of  these  parameters 

successfully,  which  are  listed  in  Table  I,  where  q^  -  Pq/)(.  this  choice, 

we  have  obtained  the  sound  velocity  c^  -  164.4  m/s  at  absolute  zero 

10 

temperature.  However,  from  the  analysis  of  the  sound  attenuation,  we 

obtained  c^  -  84.06  m/s,  which  is  very  close  to  the  experimental  value  of  76  ± 

19 

2  m/s  given  by  Wushburn  et  al,  and  the  parameters  u,  a  and  b  in  Eqs.  (2.7)  , 

/.  1  c  ^v9  21 

(2.8)  and  (2.22)  are  assumed  to  be  1.8,  1  x  10  ,  and  4  x  10  ,  as  used 

by  previous  workers  for  the  bulk  case.  With  the  use  of  the  above  parameters 
and  those  in  Table  I,  the  numerical  expressions  for  the  kinetic  coefficients 
are  given  by 


39  5 

-  1.24  x  10  T 


(4.4) 
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r5FF  ‘  l'°  *  l°43  ^ 


(4.5) 


r  ,  =  4  x  LO49  e-8-24/T 

ph-  r 


(4.6) 


Through  Eqs.  (4.4) -(4. 6)  we  see  chat  the  powers  of  the  temperatures  in  the 
kinetic  coefficients  vary  according  to  the  interactions  between  elementary 
excitations . 

Substituting  Eqs.  (4.1)-(4.4)  into  Eqs.  (  3  .  29  )  -  (  3  .  3  2 )  ,  we  obtain  the 
temperature  variation  of  the  four  coefficients  of  the  2D  second  viscosity  as 
follows : 


(T) 


4.32  x  10‘4  T" 1 f 1 


2(<7+^TZ) 


3.59  x  IQ'2  T3/2e4'12/Ti 


+  7 .19  xlO’  2  TJ/ 


3/2  4.12/T 


+  3.6] 


(4.7) 


r2(T)  - 


- 13  - 1 

2.23  x  10  u  T  1[1 


2((^)  +  3.59  x  IQ'2!  V2.4.12/T, 

- - - k - - - —  +  3.6P 


3  ^16,17  +  7.19  x10'2  T3/2e4-12/T 

4  T  T2 


(4.8) 


C3(T)  -  8.39  x  105  T'1 


(4.9) 


f4(T)  “  f1(T) 


(4.10) 


When  T  drops  to  less  than  0.3  K,  the  second  (fractional)  term  inside  the 
square  brackets  in  Eqs.  (4.7)  and  (4.8)  approaches  unity,  and  these  equations 
are  reduced  to  the  following  two  equations,  respectively: 
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r1(T)  -  1.555  x  10'3  T'1  ,  (4.11) 

f2(T)  -  2.890  x  10'12  T'1  (4.12) 


All  four  coefficients  have  a  T  3 -dependence ,  which  is  the  same  as  that  of  the 

coefficient  of  the  first  viscosity.  Here,  we  note  that  there  are  generally 

22  23 

thirteen  independent  dissipative  coefficients,  ’  but  we  only  consider  the 
above  five  coefficients.  Moreover,  we  will  not  discuss  the  pressure  and 
frequency  variations  of  the  coefficients  of  the  second  viscosity. 

We  have  already  shown  that  the  2D  thermal  conductivity  and  first 
viscosity  of  thin  liquid  helium  films33  are  all  positive.  We  can  easily 

confirm  that  the  four  coefficients  of  second  viscosity  are  all  positive.  For 

2  2-3  -2.7 

T  =«  0.05  K  we  find  that  <  C  2  5*  3  •  w^ere  -  10  and  f  ^  ^  -  10 

However,  as  the  temperature  increases  to  about  0.1  K,  the  equality  then  holds, 

2 

i.e.  ,  =«  5"  2^  3  •  Comparing  the  order  of  magnitude  for  the  four  coefficients, 

we  can  write  f  ^  f  ^  •  and  all  four  coefficients  have  larger  values  than 

that  of  the  coefficient  of  the  first  viscosity. 

In  conclusion,  the  coefficients  of  the  second  viscosity  in  two- 
dimensional  thin  helium  films  behave  like  those  of  bulk  liquid  helium  for  T  > 
0.3  K,  while  for  T  <  0.3  K,  the  four  coefficients  of  the  second  viscosity 
exhibit  a  T  3 -dependence  like  that  of  the  first  viscosity,  which  is  due  to  the 
3PP  originating  from  the  anomalous  excitation  spectrum  of  two-dimensional 
liquid  helium  at  low  momenta  and  low  temperatures. 
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Table  I.  Theoretical  parameters. 

n  (A'2) _ A/kB  (K) _ qQ  (A'1) _ Mq _ cq  (m/s> 

2.79  x  10"2  4.12  1.02  0.75 


164.4 


Figure  Capr.ion 


Figure  1.  Temperature  variation  of  the  coefficients  C ^ f 
the  second  viscosity  and  the  coefficient  rj  of  the  first 
in  thin  helium  film. 
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